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Josephson effect in a pianar graphene junction is studied by assuming that the coupling of 
a graphene sheet and two superconductors deposited on its top is described by a tunneling 
Hamiltonian. This model properly takes account of the proximity effect characteristic to a 
planar junction, and allows us to treat monolayer and bilayer cases in a parallel manner. 
Applying a quasiclassical Green's function approach to it we analyze the Josephson critical 
current I c in a short-junction limit. As a characteristic feature of the planar junction we hnd 
that J c is a concave function of temperature at the strong coupling limit while it crosses over 
to a convex function with decreasing the coupling strength. We also find different chemical- 
potential dependences of I c in the monolayer and bilayer cases. 
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Since the realization of a monolayer sheet of graphene, ^ extensive studies have been 
devoted to uncovering unusual electronic properties of this material. 2 > Josephson effect in 
graphene has been a target of intense theoretical 3-9 - 1 and experimental 10-13 ) studies during 
the last few years. The main interest is focused on how the Josephson current is affected by 
the unique band structure of graphene, i.e., the conduction and valence bands touch conically 
at K + and K- points in the Brillouin zone, and the density of states vanishes at the energy 
of the band touching point (i.e., Dirac point), which is set as e = hereafter. Titov and 
Beenakker 4 ) calculated the Josephson current through a monolayer graphene sheet on which 
two superconducting electrodes are deposited with separation L, under the assumption that 
carriers are heavily doped in the region covered by the superconductors. In the short-junction 
limit where L is much shorter than the superconducting coherence length £, they obtained the 
critical current I c at zero temperature as a function of the chemical potential [i. It is shown 
that I c is finite even at [i = and linearly increases with increasing \x. An experimental result 
consistent with this prediction has been reported. 10 ) 

We focus on another interesting aspect of the Josephson effect in graphene, stemming from 
the fact that graphene is a unique realization of an isolated ideal two-dimensional (2D) electron 
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system. Graphene in a Josephson junction acquires a 2D (planar) contact with a superconduc- 
tor, since a natural way to create a superconductor-graphene-superconductor junction is to 
deposit superconducting electrodes on top of a graphene flake. 10 ~ 13 ) This is quite contrasting 
to the case of a usual 2D electron gas imbedded in a semiconductor hetero-structure that has 
a one-dimensional (linear) contact with superconducting electrodes. However, the previous 
theoretical studies have not paid attention to the structure of such a planar junction. In the 
model used so far, 4,5,7,9 ) an energy- independent effective pair potential Aq is induced inside 
the graphene sheet over the region covered by the superconductors. This assumption reduces 
the planar junction to a conventional linear junction model. It is questionable whether the 
superconducting proximity effect in a planar junction is fully described by a conventional 
model with Aq. 

In this letter we study the stationary Josephson effect in the planar junction of graphene 
by employing a simple model in which a graphene sheet is coupled with superconductors by 
a tunneling Hamiltonian. 14 ) This model properly describes the proximity effect characteristic 
to a planar junction, which is ignored in the previous theoretical studies, allowing us to treat 
monolayer and bilayer cases in a parallel manner. We apply a quasiclassical Green's function 
approach 15 ' 16 ^ to our model under an effective mass approximation. 17 ) When fi is away from 
the Dirac point, our approach enables us to derive a general expression for the Josephson cur- 
rent, which is applicable for an arbitrary coupling strength T between the graphene sheet and 
the superconductors. Using this expression we calculate the Josephson current in the short- 
junction limit. We show that I c is a concave function of temperature T at the large-r limit, 
while it crosses over to a convex function with decreasing T. This convex T-dependence was 
not observed in the previous study 9 ) based on the model with an energy-independent effective 
pair potential, and should be regarded as a characteristic feature of the planar junction. We 
also show that the /x-dependence of I c at T = qualitatively differs in the monolayer and 
bilayer cases. We set fee = fr = I throughout this letter. 

Let us consider a clean graphene sheet on which two superconductors, Si and S2, of width 
W are deposited with separation L, where Si and S2 occupy the region of L/2 < x and 
that of x < —L/2, respectively (see Fig. 1). Note that only the top layer is in contact with 
Si and S2 in the bilayer case. Under the condition of W S> L, we regard our system as 
being translationally invariant in the y-direction. We assume that carrier doping is uniform in 
graphene and that the pair potential A(x) for the superconductors is given by A(x) = Ae 1 ^ 2 
in Si and Ae -1 ^/ 2 in S2. We use a tunneling Hamiltonian to describe the coupling of the 
graphene sheet and the superconductors. The resulting proximity effect on quasiparticles 
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Fig. 1. Planar junction geometry: Josephson junction consisting of a graphene sheet on which two 
superconductors Si and S2 of width W are deposited with separation L. The pair potential is 
assumed to be Ae 1 ^/ 2 in Si and Ae -1 ^/ 2 in S2. 



in graphene is described by a self-energy 14 ) for the thermal Green's function given below. 
We consider only quasiparticle states near the K + point because the K + and K_ points 
are degenerate. To describe quasiparticle states in graphene on the basis of a tight-binding 
model, 17,18 ) we introduce nearest-neighbor transfer integral 70. In addition, we employ nearest- 
neighbor interlayer coupling 71 and next nearest-neighbor interlayer coupling 73 in the bilayer 
case. They are estimated as 70 ~ 2.8 eV, 71 ps 0.4 eV, and 73 « 0.3 eV. 2 ) Let us introduce 
the thermal Green's function Gj(r,r';oj n ) with the Matsubara frequency co n = (2n + l)nT. 
The subscript j = 1 (j = 2) specifies the monolayer (bilayer) case. Within an effective mass 
approximation, the Green's function obeys 

(iw n f| x4 - Hj - tj) Gj{r, r'; uj n ) = f2 x4 5(r - r') (1) 

with f| x4 = diag(l, 1, —1, —1), f^ = diag(l, 1, 1, 1), Hj = diag (Hj,Hj). The 2x2 effective 
Hamiltonian Hj for low-energy quasiparticles is given by 

ft = ( 7 ) (2) 

for the monolayer case 17 ) with 7 = (\/3/2)7oa (a: lattice constant) and k± = —\d x ± d y , and 

/ -n -ak\ -$k-\ 

H2 = . (3) 

y —ak_ - /3k + -fj, J 

for the bilayer case 19 ) with a = 7 2 /7i and j3 = (y/3/2)^a characterizing the trigonal warping. 
The self-energy Ej is given by 14 ) 

7 , ~ ir ( A(X> ^ |«(W-L/2), (4) 
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where T characterizes the coupling strength between the graphene sheet and the supercon- 
ductors, X2x2 = diag(l,l), and X2x2 = diag(l,0). The matrix form of X2I2 reflects the fact 
that only the top layer is in contact with the superconductors in the bilayer case. 20 ) The off- 
diagonal elements of Sj are regarded as an energy-dependent effective pair potential, while 
the diagonal elements describe renormalization of a quasiparticle energy. If the w„-dependence 
is ignored by setting ui n = 0, our model is reduced to the conventional one. 4 ' 5 ' 9 ) 

Hereafter we restrict our attention to the moderate doping regime of 70, 71 S> p S> A. To 
introduce the quasiclassical Green's function, we perform a Fourier transformation as 

^■(p,r;w n ) = J tfse-^'&j (r + |,r - |;w n ) , (5) 

where p = (p x ,Py)- In the monolayer case of j = 1, the Green's function obeys 

(iw n f| x4 - U\- Si) Gi{p,r-u n ) = ff x4 , (6) 

where the 4x4 Hamiltonian iii is given by iii = diag(i?i + hi, Hi + hi) with 

Hi=\ -» T^-'vyA, (7) 

l(Px + iPy) ~P J 

2jfe_ \ 

p+ y 

The 2x2 matrix Hi is diagonalized as u P Hiu p = diag (jp — p, —jp — p) in terms of 

where p = \p\ and </> p = aig{p x + \p y }. As long as // 3> A, the subband with the energy disper- 
sion — 7j> — is irrelevant in the superconducting proximity effect. Therefore, we are allowed 
to consider only the relevant subband with the energy dispersion 7^5 — p. In accordance with 
this observation, we transform Gi as Qi(p, r; ui n ) = U p Gi{p, r; u n )U p with U p = diag (u p , u p ) 
and retain only the (1, 1)-, (1,3)-, (3, 1)-, and (3, 3)-elements. 20 ) Accordingly, we define Gi as 

Gi=( ^ ^1.3 V (10) 

The Green's function Gi approximately satisfies 

iw n Tf X 2 - (lP- M + ^vi(p) • V^j - Si Gi(p,r;u n ) = t$ x2 , (11) 
where r| x2 = diag(l, — 1), r ^ = diag(l,l), the velocity «i(p) is given by «i(p) = 
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7(cos (j) p , sin <pp) , and 

_ -if I u n A(x) 
V A2 + < \ A{xf -u> n 
By repeating the argument similar to this we can show that the Green's function for the 
bilayer case satisfies 



9(\x\-L/2). 



(12) 



iw n r 2x 2 



G 2 (p,r;u; n ) = r 2 ° x2 , 



(13) 



where e p = [(ap 2 ) 2 + 2af3p 3 cos(30 p ) + (/3p) 2 ] 1//2 , and £ 2 = Ei/2. The x- and y-components 
of vzip) are given by 

V2x(p) = 2apcos(9 p - 4> p ) + (3 cos 9 P , (14) 

V2y(p) = 2ap$m(9p - cp p ) - f3s'm9 p , (15) 

where 9 P = arg{ap 2 e l2 ^ p + (3pe~ l< t )p }. Note that S 2 is smaller by a factor of two than Si 
reflecting the fact that only the top layer is in contact with the superconductors. 
We define the quasiclassical Green's function Gj(n,r;ui n ) with n = p/p as 21 ' 22 ) 

(16) 



7T 



Gj(n,r;u n ) = - d^pGj(p,r;u} n ), 



where a diverging contribution must be subtracted, and £ p = jp — fi for j = 1 and £ p = e p — (i 
for j = 2. Applying a standard procedure 22 ) to eqs. (11) and (13), we can show that the 
quasiclassical Green's function satisfies 



' luJ n [r| x2 , Gj] + iv Fj (n) ■ VGj - Gj] = 0, 



(17) 



where vpj(n) represents the Fermi velocity in the momentum direction denoted by n. We 
express the elements of Gj(n,r;u n ) as [Gj]i,i = -[Gj] 2 ,2 = 9j, [Gj]i,2 = fj, and [Gj] 2 ,i = /j. 
Equation (17) yields 



v Fj (n) ■ V 9j = Q(x,u> n ) A( X yfj - A(x)f) 



2 [1 + (j(x,uj n )} LUnfj + v F j(n) ■ Vfj = 2( j (x,u n )A(x)g j , 
2 [1 + (j(x,uj n )} u) n f] - v Fj (n) ■ V/J = 2Q(x, uj n )A(x)*gj, 



where 



\x\-L/2) 



(18) 
(19) 
(20) 

(21) 



with Ti = T and T 2 = T/2. Note that although we are considering a Josephson junction of 
graphene, only the Fermi velocity reflects a feature of graphene in eqs. (18)-(20). We point 
out that these equations are also applicable to a conventional planar junction of a 2D electron 
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system. 23 ^ 

We derive a general expression for the dc Josephson current Ij(<p) in terms of the qua- 
siclassical Green's function, where j = 1 and 2 correspond to the monolayer and bilayer 
cases, respectively. We ignore the trigonal warping effect by setting ft = 0. With this sim- 
plification, the Fermi velocity is expressed as v^j = vpj (cos (ft, sin (ft) with vpi = 7 and 



vp2 = 27yV/7i- The density of states per spin at the Fermi level is given by -/Vi(O) = fi/ (ttj 2 ) 
and iV2(0) = 7!/(27r7 2 ), in which the two-fold valley degeneracy is included. The Josephson 
current is expressed as 15 ' 16 ) 

Ij(tp) = 2^(0) W 1 2 ^v Fjx TJ2lm{ gj (n,x;cu n )} . (22) 

Note that in the r-dependence of gj, only x is relevant since VF » L is assumed. Because 
Ij((/j) is independent of x for \x\ < L/2, we evaluate it at x = L/2. Parameterizing x as 
x = (vpj cos(ft)t, we solve eqs. (18)-(20) with gj + fjfj = 1 under the boundary condition of 

lim gj(n,t;uj n ) = ^, (23) 

Ae ±i ^ 2 

lim fj(n,t;u n ) = ~- , (24) 

t->±oo 



where w = [1 + (j(x, co n )] co n , A = (j(x,co n )A, and 17 = \/A 2 + w 2 . Note that x = L/2 

corresponds to t+ = L/{2v-pj cos (ft). A straightforward calculation yields 

Co A 2 (sinh k cos % + icosh Ksin %) , 

g j (n,t+;u; n ) = 7 - + —r Z * , 2 . , ^ r (25) 

s ' 17 ( Q cosh k + Co sinh k J cos ^ + iQ IQ, sinh k + oj cosh k J sin |- 

with /c = w n L/ (vpj cos </>). We obtain 

i iM = 2^,(0)^ /» ^*rE 7= \ ^ ^ — ■ (26) 

J-f n Un [Q? + w 2 J cosh 2k + 2wOsinh2K + A 2 cos 99 

Using this general expression one can numerically calculate the Josephson current in the 
planar junction for arbitrary parameters. In the strong-coupling limit of Tj S> Ao with Ao 
being the magnitude of the pair potential at T = 0, eq. (26) is reduced to the ordinary 
expression for the Josephson current through a 2D electron gas of finite area placed between 
two superconductors. 15 ' 24 ^ 

Below we focus on the short-junction limit of L <C £, where £ = VFj/(27rAo) is the super- 
conducting coherence length, and study the behavior of the critical current for an arbitrary 
Tj. We can approximate as cosh 2k « 1 and sinh 2k rj in the short-junction limit, and obtain 

„ r 2 A 2 sinw 

Ij (99) = eAfjT Y 3 r , (27) 

tT u\ (A 2 + a; 2 + 2r, v^T^ 2 + T 2 ) + T 2 A 2 cos 2 I 
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where Mj = 2vpjNj{fS)W represents the number of conducting channels. Before considering 
the critical current, let us observe the behavior of Ij(ip) in the strong and weak coupling limits. 
Firstly we consider the strong coupling limit of S> Aq > A. In this case we can ignore all 



terms in the parentheses except for V 2 , and obtain 



= eftfj A sin ^ tanh ^ 



A cos ^ 



Ij((p) = eMjTj sin £ tanh ( J orTi 2 ) . (29) 



2T " (28) 
This is identical to the result derived by Kulik and OmeFyanchek (KO). 16 ) We next consider 
the weak coupling limit of Ao » Tj. In the low-temperature regime of A » Fj » T, we can 
ignore all terms in the parentheses except for A 2 , and obtain 23 ) 

^tanhf^ 2 
2 \ 2T 

When T c > T, there holds T » Tj, A which enables us to retain only ui^ in the denominator. 
We thus obtain 

r 2 A 2 

I j {p>) = eN j ^smp>. (30) 

This T-dependence is qualitatively different from the KO result which yields I(tp) = 

eMA 2 (4T)- 1 simp near T c . 

To observe the T-dependence of the critical current 1^ = max^ {Ij((p)}, we numerically 

calculate 1^ for several values of r = Tj/Ao on the basis of eq. (27). The T-dependence of 

A is determined by the gap equation 1 = Amt Jq D de tanh( y/<: ^+ A -)/ Ve 2 + A 2 , where Ai nt is 

the dimensionless interaction constant, and the Debye energy is chosen as €d/Ao = 200. The 

critical current 1^ normalized by 1^ = eNjAo is shown in Fig. 2 as a function of T/T c for 

r = 0.2, 1.0, 5.0, and r — > oo at which the KO result is reproduced. Note that 1^ /I^ does 

(i) 

not depend on j. We observe that Ic is a concave function of T at r — > oo, while it crosses 
over to a convex function with decreasing r. Such a convex T-dependence was not observed 
in the previous study 9 ) based on an energy-independent effective pair potential model. The 
coupling strength Tj crucially affects the T-dependence of the critical current. 

Let us finally evaluate the critical current at T = 0. As an intermediate step, we derive 
a convenient analytic formula by approximating the w n -dependence of the denominator of 



eq. (27) as 2Tj \J A 2 + uj 2 ~ 2TjA. This is justified either in the low-temperature limit of 
T c S> T for an arbitrary Tj, or in the strong and weak coupling limits for an arbitrary T. 
Performing the summation over uj n in terms of a contour integral, we obtain 



T 2 A 2 sin ip 

Ij(<p) = eAfj 



^E + {p)E„(p) 



tanh ( s+frMM?) ) tanh ( 



E+(<p)+E-{<py 

AT 



E+(<p) - E.(<p) E + (<p) + E.(<p) 



(31) 
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Fig. 2. Temperature dependence of the normalized critical current 1^/1^ for r = 0.2, 1.0, 5.0, and 
r — > oo, where r = Fj/Aq. 



where E±(ip) = [(A + Tj) 2 ± cos(( / 9/2)] 1 / 2 . One can verify that eq. (31) reproduces the 
correct asymptotic behaviors, eqs. (28)-(30), in the corresponding limits. Note that eq. (31) 
is maximized at (p = ir (mod 2tt) in the limit of T — > 0. The critical current is determined as 
= eA/" j r j A /(A + Tj), which yields 



2fiW TA 

if =e 

for the monolayer, and 



/(1)=C W 1*0 (32) 
7T 7 A + r v ; 



/(2) = g w /tT i a q 



vt7 V n A + £ 



(33) 



for the bilayer cases. We see that 1^ is proportional to /i, which is consistent with the result 

(2) 

reported in ref. 4, while Ic is proportional to y/Ji. Roughly speaking, the critical current in 
the bilayer case is greater than that in the monolayer case by a factor of y^fi/Ji- 

In summary we have proposed a model for a planar Josephson junction of graphene, and 
derived a general expression for the Josephson current at moderate doping in the quasiclassical 
Green's function approach. Much emphasis has been on the behavior of the Josephson current 
in the short-junction limit in monolayer and bilayer graphene junctions. It was demonstrated 
that the coupling strength crucially affects the temperature dependence of the critical current 
in an unexpected manner. This should be regarded as a characteristic feature of the planar 
junction. We have also shown that the chemical-potential dependence of the critical current 
qualitatively differs in the monolayer and bilayer cases. Finally we point out that our argument 
can be extended to a multilayer case. 20 ) Such an extension will be reported elsewhere. 
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